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Abstract 

A recently introduced systematic approach to derivations of the macroscopic dy- 
namics from the underlying microscopic equations of motions in the short-memory 
approximation [Gorban et al, Phys. Rev. E 63, 066124 (2001)] is presented in de- 
tail. The essence of this method is a consistent implementation of Ehrenfest's idea of 
coarse-graining, realized via a matched expansion of both the microscopic and the 
macroscopic motions. Applications of this method to a derivation of the nonlinear 
Vlasov-Fokker-Planck equation, diffusion equation and hydrodynamic equations of 
the fluid with a long-range mean field interaction are presented in full detail. The 
advantage of the method is illustrated by the computation of the post-Navier-Stokes 
approximation of the hydrodynamics which is shown to be stable unlike the Burnett 
hydrodynamics. 
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1 Introduction 



The question of how irreversibility can be derived from reversible dynamics is 
one of the classical problems in physics. The first solution has been suggested 
by Boltzmann [1], and it provoked much discussion at that time. An alternative 
approach has been given by Ehrenfest [2] who coined the notion of coarse- 
graining. 

The impact of Ehrenfest 's ideas on the long-standing discussions of the foun- 
dations of the nonequilibrium thermodynamics is enormous (see, e. g. [3,4]). 
In a recent paper [5] we have given a novel formalization of Ehrenfest's ap- 
proach. The main focus of Ref. [5] was the mathematical consistency of the 
formalization, whereas applications were only briefly indicated. The goal of 
the present paper is to give a detailed description of the method, focusing on 
how to apply it to various typical examples. 

The starting point of our construction are microscopic equations of motion. A 
traditional example of the microscopic description is the Liouville equation for 
classical particles. However, we need to stress that the distinction between "mi- 
cro" and "macro" is always context dependent. For example, Vlasov's equation 
describes the dynamics of the one-particle distribution function. In one state- 
ment of the problem, this is a microscopic dynamics in comparison to the 
evolution of hydrodynamic moments of the distribution function. In a differ- 
ent setting, this equation itself is a result of reducing the description from the 
microscopic Liouville equation. 

The problem of reducing the description includes a definition of the micro- 
scopic dynamics, and of the macroscopic variables of interest, for which equa- 
tions of the reduced description must be found. The next step is the construc- 
tion of the initial approximation. This is the well known quasi-equilibrium ap- 
proximation, which is the solution to the variational problem, S — > max, where 
S in the entropy, under given constraints. This solution assumes that the mi- 
croscopic distribution functions depend on time only through their dependence 
on the macroscopic variables. Direct substitution of the quasi-equilibrium dis- 
tribution function into the microscopic equation of motion gives the initial 
approximation to the macroscopic dynamics. All further corrections can be 
obtained from a more precise approximation of the microscopic as well as 
of the macroscopic trajectories within a given time interval r which is the 
parameter of our method. 
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The method described here has several clear advantages: 

(i) It allows to derive complicated macroscopic equations, instead of writing 
them ad hoc. This fact is especially significant for the description of complex 
fluids. The method gives explicit expressions for relevant variables with one 
unknown parameter (r). This parameter can be obtained from the experimen- 
tal data. 

(ii) Another advantage of the method is its simplicity. For example, in the 
case where the microscopic dynamics is given by the Boltzmann equation, the 
approach avoids evaluation of Boltzmann collision integral. 

(iii) The most significant advantage of this formalization is that it is appli- 
cable to nonlinear systems. Usually, in the classical approaches to reduced 
description, the microscopic equation of motion is linear. In that case, one 
can formally write the evolution operator in the exponential form. Obviously, 
this does not work for nonlinear systems, such as, for example, systems with 
mean field interactions. The method which we are presenting here is based on 
mapping the expanded microscopic trajectory into the consistently expanded 
macroscopic trajectory. This does not require linearity. Moreover, the order- 
by-order recurrent construction can be, in principle, enhanced by restoring to 
other types of approximations, like Pade approximation, for example, but we 
do not consider these options here. 

In the present paper we discuss in detail applications of the method [5] to 
derivations of macroscopic equations in various cases, with and without mean 
field interaction potentials, for various choices of macroscopic variables, and 
demonstrate how computations are performed in the higher orders of the ex- 
pansion. The structure of the paper is as follows: In section 2, for the sake of 
completeness, we describe briefly the formalization of Ehrenfest's approach [5]. 
We stress the role of the quasi-equilibrium approximation as the starting point 
for the constructions to follow. We derive explicit expressions for the correction 
to the quasi-equilibrium dynamics, and conclude this section with the entropy 
production formula and its discussion. In section 3, we begin the discussion of 
applications. The first example is the derivation of the Fokker-Planck equation 
from the Liouville equation. In section 4 we use the present formalism in or- 
der to derive hydrodynamic equations. Zeroth approximation of the scheme is 
the Euler equations of the compressible nonviscous fluid. The first approxima- 
tion leads to the system of Navier-Stokes equations. Moreover, the approach 
allows to obtain the next correction, so-called post-Navier-Stokes equations. 
The latter example is of particular interest. Indeed, it is well known that post- 
Navier-Stokes equations as derived from the Boltzmann kinetic equation by 
the Chapman-Enskog method (Burnett and super-Burnett hydrodynamics) 
suffer from unphysical instability already in the linear approximation [6]. We 
demonstrate it by the explicit computation that the linearized higher-order 
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hydrodynamic equations derived within our method are free from this draw- 
back. In section 5, we derive macroscopic equations in the case of the nonlinear 
microscopic dynamics (hydrodynamic equations from the Vlasov kinetic equa- 
tion) . Significance of this example is that methods based on the projection 
operator approach [7,8,9], are inapplicable to nonlinear systems. We show in 
detail how this problem is solved on the basis of our method. 



2 General construction 



Let us consider a microscopic dynamics given by an equation, 

f = Af), (1) 



where f(x,t) is a distribution function over the phase space x at time t, and 
where operator J(f) may be linear or nonlinear. We consider linear macro- 
scopic variables M k = //&(/), where operator /i k maps / into M k . The prob- 
lem is to obtain closed macroscopic equations of motion, M k = <f> k (M). This 
is achieved in two steps: First, we construct an initial approximation to the 
macroscopic dynamics and, second, this approximation is further corrected on 
the basis of the coarse-gaining. 

The initial approximation is the quasi-equilibrium approximation, and it is 
based on the entropy maximum principle under fixed constraints [10,11]: 

S(f) - max, = M, (2) 



where S is the entropy functional, which is assumed to be strictly concave, 
and M is the set of the macroscopic variables {M}, and \x is the set of the 
corresponding operators. If the solution to the problem (2) exists, it is unique 
thanks to the concavity of the entropy functionals. Solution to equation (2) 
is called the quasi-equilibrium state, and it will be denoted as f*(M). The 
classical example is the local equilibrium of the ideal gas: / is the one-body 
distribution function, S is the Boltzmann entropy, /i are five linear operators, 
M/) = /{!> v ' v2 }fdv, with v the particle's velocity; the corresponding f*(M) 
is called the local Maxwell distribution function. 

If the microscopic dynamics is given by equation (1), then the quasi-equilibrium 
dynamics of the variables M reads: 

M k = ^k(J(r(M)) = <f>l (3) 
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The quasi-equilibrium approximation has important property, it conserves the 
type of the dynamics: If the entropy monotonically increases (or not decreases) 
due to equation (1), then the same is true for the quasi-equilibrium entropy 
S*(M) = S(f*(M)), due to the quasi-equilibrium dynamics (3). That is, if 

* = M£> y = «|!) J(/) > , 



then 



V M k 

^dM k k 



dS 



E^W*(M)))>0. 



(4) 



Summation in k always implies summation or integration over the set of labels 
of the macroscopic variables. 

Conservation of the type of dynamics by the quasi-equilibrium approximation 
is a simple yet a general and useful fact. If the entropy S is an integral of motion 
of equation (1) then S*(M) is the integral of motion for the quasi-equilibrium 
equation (3). Consequently, if we start with a system which conserves the 
entropy (for example, with the Liouville equation) then we end up with the 
quasi-equilibrium system which conserves the quasi-equilibrium entropy. For 
instance, if M is the one-body distribution function, and (1) is the (reversible) 
Liouville equation, then (3) is the Vlasov equation which is reversible, too. On 
the other hand, if the entropy was monotonically increasing on solutions to 
equation (1), then the quasi-equilibrium entropy also increases monotonically 
on solutions to the quasi-equilibrium dynamic equations (3). For instance, if 
equation (1) is the Boltzmann equation for the one-body distribution function, 
and M is a finite set of moments (chosen in such a way that the solution to 
the problem (2) exists), then (3) are closed moment equations for M which 
increase the quasi-equilibrium entropy (this is the essence of a well known 
generalization of Grad's moment method). 



2. 1 Enhancement of quasi-equilibrium approximations for entropy- conserving 
dynamics 



The goal of the present section is to describe the simplest analytic implemen- 
tation, the microscopic motion with periodic coarse-graining. The notion of 
coarse-graining was introduced by P. and T. Ehrenfest's in their seminal work 
[2]: The phase space is partitioned into cells, the coarse-grained variables are 
the amounts of the phase density inside the cells. Dynamics is described by 
the two processes, by the Liouville equation for /, and by periodic coarse- 
graining, replacement of f(x) in each cell by its average value in this cell. The 
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coarse-graining operation means forgetting the microscopic details, or of the 
history. 

From the perspective of general quasi-equilibrium approximations, periodic 
coarse-graining amounts to the return of the true microscopic trajectory on the 
quasi-equilibrium manifold with the preservation of the macroscopic variables. 
The motion starts at the quasi-equilibrium state /*. Then the true solution 
fi(t) of the microscopic equation (1) with the initial condition fi(0) = f* is 
coarse-grained at a fixed time t — r, solution fi(r) is replaced by the quasi- 
equilibrium function f* +1 = /*(/i(/j(r))). This process is sketched in Fig. (1). 

From the features of the quasi-equilibrium approximation it follows that for 
the motion with periodic coarse-graining, the inequality is valid, 

sun < s(f* +1 ), (5) 



the equality occurs if and only if the quasi-equilibrium is the invariant manifold 
of the dynamic system (1). Whenever the quasi-equilibrium is not the solution 
to equation (1), the strict inequality in (5) demonstrates the entropy increase. 

In other words, let us assume that the trajectory begins at the quasi-equilibrium 
manifold, then it takes off from this manifold according to the microscopic 
evolution equations. Then, after some time r, the trajectory is coarse-grained, 
that is the, state is brought back on the quasi-equilibrium manifold keeping 
the values of the macroscopic variables. The irreversibility is born in the lat- 
ter process, and this construction clearly rules out quasi-equilibrium manifolds 
which are invariant with respect to the microscopic dynamics, as candidates 
for a coarse-graining. The coarse-graining indicates the way to derive equations 
for macroscopic variables from the condition that the macroscopic trajectory, 
M{t), which governs the motion of the quasi-equilibrium states, f*(M(t)), 
should match precisely the same points on the quasi-equilibrium manifold, 
f*(M(t + r)), and this matching should be independent of both the initial 
time, t, and the initial condition M(t). The problem is then how to derive the 
continuous time macroscopic dynamics which would be consistent with this 
picture. The simplest realization suggested in the Ref. [5] is based on using 
an expansion of both the microscopic and the macroscopic trajectories. Here 
we present this construction to the third order accuracy, in a general form, 
whereas only the second-order accurate construction has been discussed in [5] . 

Let us write down the solution to the microscopic equation (1), and approxi- 
mate this solution by the polynomial of third oder in r. Introducing notation, 
J* = J(f*(M(t))), we write, 

2Q7* 3 /f)J*f)J* fp j* \ 
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Evaluation of the macroscopic variables on the function (6) gives 



M k (t + r) = M k + r<t>l + -fi k [j^ J *j 

(7 ) 

where (p k = fi k (J*) is the quasi-equilibrium macroscopic vector field (the right 
hand side of equation (3)), and all the functions and derivatives are taken in 
the quasi-equilibrium state at time t. 

We shall now establish the macroscopic dynamic by matching the macro- 
scopic and the microscopic dynamics. Specifically, the macroscopic dynamic 
equations (3) with the right-hand side not yet defined, give the following third- 
order result: 



M k {t + T)=M k + T<p k + -Y J -^<t> j 

3 3 

+ 4E( /^-M- + l^l^l + o(r% (8) 
3! Y \dMiMj J dMi dMj 3 J 1 ; 1 ; 



Expanding functions (f) k into the series (f) k = R k + rR k + T 2 R k + 
(R k °^ = 0*), and requiring that the microscopic and the macroscopic dynamics 
coincide to the order of r 3 , we obtain the sequence of corrections for the right- 
hand side of the equation for the macroscopic variables. Zeroth order is the 
quasi-equilibrium approximation to the macroscopic dynamics. The first-order 
correction gives: 



The next, second-order correction has the following explicit form: 



(2 ) _ 1 f (dJ* d.P \ (d 2 J* \\ i f d<pi d<p* \ 

- gj E [e^W) - a E [swA + . 

Further corrections are found by the same token. Equations (9)-(10) give 
explicit closed expressions for corrections to the quasi-equilibrium dynamics 
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to the order of accuracy specified above. They are used below in various specific 
examples. 

2.2 Entropy production 

The most important consequence of the above construction is that the result- 
ing continuous time macroscopic equations retain the dissipation property of 
the discrete time coarse-graining (5) on each order of approximation n > 1. 
Let us first consider the entropy production formula for the first-order ap- 
proximation. In order to shorten notations, it is convenient to introduce the 
quasi-equilibrium projection operator, 

^=e|£ms). (id 

It has been demonstrated in [5] that the entropy production, 




equals 



(1-P*)J*. (12) 

/• 

Equation (12) is nonnegative definite due to concavity of the entropy. Entropy 
production (12) is equal to zero only if the quasi-equilibrium approximation 
is the true solution to the microscopic dynamics, that is, if (1 — P*)J* = 0. 
While quasi-equilibrium approximations which solve the Liouville equation are 
uninteresting objects (except, of course, for the equilibrium itself), vanishing of 
the entropy production in this case is a simple test of consistency of the theory. 
Note that the entropy production (12) is proportional to r. Note also that 
projection operator does not appear in our consideration a priory, rather, it is 
the result of exploring the coarse-graining condition in the previous section. 

Though equation (12) looks very natural, its existence is rather subtle. In- 
deed, equation (12) is a difference of the two terms, J2k Hk(J*dJ* / 'df) (contri- 
bution of the second-order approximation to the microscopic trajectory), and 
J2ik dR^ /dMi (contribution of the derivative of the quasi-equilibrium vec- 
tor field). Each of these expressions separately gives a positive contribution to 
the entropy production, and equation (12) is the difference of the two positive 
definite expressions. In the higher order approximations, these subtractions are 



% = -i(i-n/-f! 
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more involved, and explicit demonstration of the entropy production formulae 
becomes a formidable task. Yet, it is possible to demonstrate the increase- 
in-entropy without explicit computation, though at a price of smallness of r. 
Indeed, let us denote S* n ^ the time derivative of the entropy on the nth order 
approximation. Then 

t+T 

J S{ n) (s)ds = S*(t + r) - S*(t) + 0(r n+1 ), 
t 

where S*(t + r) and S*(t) are true values of the entropy at the adjacent states 
of the if-curve. The difference 5S = S*(t + r) - S*(t) is strictly positive for 
any fixed r, and, by equation (12), 5S ~ r 2 for small r. Therefore, if r is small 
enough, the right hand side in the above expression is positive, and 

rS* {n) (9 {n) ) > 0, 

where t < 0( n ) < t + r. Finally, since Sf n Jt) = SLJs) + 0(r n ) for any s on the 

segment [t,t + r], we can replace (#(«)) in the latter inequality by Sf n Jt). 
The sense of this consideration is as follows: Since the entropy production 
formula (12) is valid in the leading order of the construction, the entropy 
production will not collapse in the higher orders at least if the coarse-graining 
time is small enough. More refined estimations can be obtained only from the 
explicit analysis of the higher-order corrections. 

2. 3 Relation to the work of Lewis 

Among various realizations of the coarse-graining procedures, the work of 
Lewis [12] appears to be most close to our approach. It is therefore pertinent 
to discuss the differences. Both methods are based on the coarse-graining 
condition, 

M k (t + r)=fi k (Trf*(M(t))), (13) 

where T T is the formal solution operator of the microscopic dynamics. Above, 
we applied a consistent expansion of both, the left hand side and the right hand 
side of the coarse-graining condition (13), in terms of the coarse-graining time 
r. In the work of Lewis [12], it was suggested, as a general way to exploring 
the condition (13), to write the first-order equation for M in the form of the 
differential pursuit, 

Mk (t) + r^p- » u fc (T T f*(M(t))) . (14) 
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In other words, in the work of Lewis [12], the expansion to the first order was 
considered on the left (macroscopic) side of equation (13), whereas the right 
hand side containing the microscopic trajectory T T f*(M(t)) was not treated 
on the same footing. Clearly, expansion of the right hand side to first order 
in t is the only equation which is common in both approaches, and this is 
the quasi-equilibrium dynamics. However, the difference occurs already in the 
next, second-order term (see Ref. [5] for details). Namely, the expansion to 
the second order of the right hand side of Lewis' equation (14) results in a 
dissipative equation (in the case of the Liouville equation, for example) which 
remains dissipative even if the quasi-equilibrium approximation is the exact 
solution to the microscopic dynamics, that is, when microscopic trajectories 
once started on the quasi-equilibrium manifold belong to it in all the later 
times, and thus no dissipation can be born by any coarse-graining. 

On the other hand, our approach assumes a certain smoothness of trajectories 
so that application of the low-order expansion bears physical significance. For 
example, while using lower-order truncations it is not possible to derive the 
Boltzmann equation because in that case the relevant quasi-equilibrium mani- 
fold (iV-body distribution function is proportional to the product of one-body 
distributions, or uncorrelated states, see next section) is almost invariant dur- 
ing the long time (of the order of the mean free flight of particles), while the 
trajectory steeply leaves this manifold during the short-time pair collision. It 
is clear that in such a case lower-order expansions of the microscopic trajec- 
tory do not lead to useful results. It has been clearly stated by Lewis [12], that 
the exploration of the condition (13) depends on the physical situation, and 
how one makes approximations. In fact, derivation of the Boltzmann equation 
given by Lewis on the basis of the condition (13) does not follow the differential 
pursuit approximation: As is well known, the expansion in terms of particle's 
density of the solution to the BBGKY hierarchy is singular, and begins with 
the linear in time term. Assuming the quasi-equilibrium approximation for 
the A r -body distribution function under fixed one-body distribution function, 
and that collisions are well localized in space and time, one gets on the right 
hand side of equation (13), 

f(t + T) = f(t)+nTJ B (f(t)) + o(n), 

where n is particle's density, / is the one-particle distribution function, and 
Jb is the Boltzmann's collision integral. Next, using the mean-value theorem 
on the left hand side of the equation (13), the Boltzmann equation is derived 
(see also a recent elegant renormalizat ion-group argument for this derivation 
[13]). 

We stress that our approach of matched expansion for exploring the coarse- 
graining condition (13) is, in fact, the exact (formal) statement that the un- 
known macroscopic dynamics which causes the shift of M k on the left hand 
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side of equation (13) can be reconstructed order-by-order to any degree of ac- 
curacy, whereas the low-order truncations may be useful for certain physical 
situations. A thorough study of the cases beyond the lower-order truncations 
is of great importance which is left for future work. 



3 Vlasov-Fokker-Planck kinetic equation 



In this section we derive kinetic equations based on the approach formulated 
above. Here microscopic dynamics is given by the N— particle Liouville equa- 
tion. Macroscopic variable is the one-particle distribution function. The so- 
lution to the variational problem (2) is the approximation of the absence of 
correlations. In this case the quasi-equilibrium iV-particle distribution func- 
tion is proportional to the product of the one-particle distribution functions. 
On the basis of this quasi-equilibrium we obtain the Vlasov equation, as the 
zeroth approximation, and the Fokker-Planck equation, as the next correction. 

The dynamics of the ensemble of N classical point particles interacting by a 
pair potential is given by the Liouville equation, 



Lw N = {//. w N } 



dt 



/ <9$(|rj - rj\) Own | dU(rj) dw N pj dw N \ 
" '• <)r, dpi dvi dpi in i)r, ) ' 



ij \ 

where wn is the N— particle distribution function, 



w N = w N (pi, ...pat, ri, ...r N , t). 

It is normalized to one (/ w N dp 1 ...dr N = 1). L is the Liouville operator (Pois- 
son bracket), H is classical Hamilton function, $(|rj— 1\,-|) describes interaction 
between particles with indices % and j, U(ri) is an external field. Here, and 
in every case later, summation is assumed over the total number of particles. 
The indices i and j take values from 1 to N independently. 

The macroscopic variable is the one-particle distribution function f(p,r,t). 
The latter satisfies the normalization condition / f(p,r,t)drdp = N. The 
mapping of the microscopic variables into the space of macroscopic variables 
is given by the following operator: 



H = J2 J $( x - Vi)dyi...dy N , 



(16) 
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where x denotes the point (p, r) in the phase space. The solution to the vari- 
ational problem (2) has the form: 



w* N ( Xl , x N , t) = N~ N J] f( Xi , t), (17) 



Before proceeding further, a remark on the choice of the one-particle distri- 
bution function as the macroscopic variable is in order. The total energy of 
the system under consideration can be expressed as a linear functional of the 
two-particle distribution function rather than in terms of the one-particle dis- 
tribution. For that reason, when only the one-particle distribution is chosen 
for the macroscopic variable, the coarse-graining procedure must be supple- 
mented by a termostatting procedure in order to keep the energy balance in 
the system intact. While it is possible to take care of the energy conserva- 
tion through introducing extra terms into the Liouville equation describing 
interactions with the thermostat, we will implement thermostatting after the 
coarse-graining. Now we proceed with executing the formalism developed in 
the previous sections. 

The quasi-equilibrium (conservative) dynamics in the space of the macroscopic 
variables is given by the action of the operator (16) on the microscopic equation 
of motion (15). 



df(x,t) 
dt 



J2j 5 ( x ~ yi)dyi---dy N 



| 9U(r] k ) dw^ _ p k dw* N ^ <9$(|% - r)i\) dw* N 
2-^ fi^. pit. ™ fi^. 



k 



dr] k d£ k m drj k f^ k dr) k d£ k 



In order to avoid a confusion, we used variables rj, £, and y under integrals 
instead of r, p, and x, respectively. Thus, we obtain: 

df(x,t) + p df(x,t) _ dUjv) df(x,t) + F , t) df(x,t) = ^ 
dt m dr dr dp dp 

where 

F(rtt)= jmtp\l fMdy (19) 

is the mean field interaction force. Equation (18) has been first derived by 
Vlasov, and is usually applied to description of plasma without collisions 
[14,15]. 
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Now we are going to obtain dissipative correction to the Vlasov equation (18) 
based on the approach developed above. For the sake of simplicity we omit 
external field term. Let us begin with the second term in the expression (9). 



d(j)* 
dM 



P\ 2 d_fdf\ 
my dr \dr J 



where 



(20) 



m J or \ or / 



(21) 



First term in equation (9) is proportional to: 



J* = L 2 w* 



N 



y-_Pij9_ + y- 0$(|r,-rj|) d 
/ — m <)r, ,4^, <)v, dpi 



^_Pk_d_ + ^d$(\r k -r l \) d 



m dru 

k K k,ly±l 



xHf(x s ,t)N- N , 



dr h 



dp k 



(22) 



where indices i, j, k, and / enumerate particles in the system. 
Removing the brackets in (22), we obtain, 



v P*Pfc d 2 2V p ^$(|r fc -r,|) d 2 v 1 d^dn-Tjl) d 



E 

ijkl 



m 2 dridr k fj^ k m dr k 
^(Iri-rj-D^drfc-r,!) d 2 



dridp k ^ m dri 



dpi 



dr k 



dr k 



dpidp k 



Uf(x s ,t)N- N , 



(23) 



where j ^ i and k ^ I. 



Let us now consider the first term in (23) under the action of the operator 
(16). The result is not equal to zero only for % — k. In this case we have: 
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\mj or \or J 

Note that this term is cancelled by the same term in equation (20). The second 
and the third term in equation (23) do not contribute to the final result by 
the same reason. 

Let us consider the last term in (23), whose contribution is nontrivial. It is: 

^(Iri-rj-D^drfc-r,!) 



2 / s ( x -y P )J2 
d 2 



Hf(y s ,t)N- N d yi ..dy N , (24) 



where k ^ I and i ^ j. 

i) For j = I, expression (24) gives: 



dp 2 J \ dv J J y ' N N 

In the thermodynamic limit (N — > oo) the expression (25) becomes equal to 
the following: 



ii) For j ^ I, expression (24) gives: 



*f(f^\r-rV md ,YN^^. (2?) 



dp 2 \J dr J v ' J N N ^ 

Combining the expressions (26) and (27) with the last term in the equation 
(20), we obtain: 



(-^)g(/^'H'- 

In the thermodynamic limit the term 1/N 2 can be neglected. 
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Thus, we obtain the following macroscopic equation of motion: 



= D 




f(y,t)dy 



(28) 



where D is the diffusion tensor, 




d<5>(\r-ri\) 
dr 



f(y,t)dy) {29) 




Note the clear Green-Kubo-type structure of the latter expression. The ab- 
sence of the time integration reveals the short-memory nature of the construc- 
tion. 

At this point of the derivation, one can notice that the obtained collision in- 
tegral in equation (28) does not conserve the total energy. As we have argued 
above, this fact has the following explanation: Since the interaction between 
particles is pair-wise, the total energy is a functional of the two-particles dis- 
tribution function. However, we have restricted ourselves to the equation for 
the one-particle distribution function function. At the beginning of our pro- 
cedure, the trajectories belong to the quasi-equilibrium manifold, later the 
correlations start growing, and lead to the decrease of the potential energy, 
and to the increase of the kinetic energy. At every instance of the microscopic 
motion, the total energy is conserved. At some moment (defined by the step of 
coarse-graining) the system is returned back onto the quasi-equilibrium mani- 
fold with the conservation of value of all the variables which can be expressed 
as functionals of the one-particle distribution function (including the kinetic 
energy). As a result, the total energy is not conserved. There are two ways to 
solve this problem. The first is to choose the two-particle distribution function 
as the macroscopic variable. This route is very complicated, and in order to 
circumvent it and to stay on the level of the one-particle distribution func- 
tion, we shall subtract the spurious contribution by a regularization of the 
macroscopic vector field. 

Specifically, let us write the dissipative contribution to equation (28) in the 
gradient form: 




9j_ 



(30) 



here 
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This form automatically takes into account the conservation of particle's den- 
sity. In order to satisfy the conservation of the momentum and of the energy, 
we shall introduce a subspace E, and require that j belong to E. Subspace E is 
a set of functions {</?} for which the momentum and the energy are conserved. 
This means that the moments, 




are equal to zero (M = 0, M a = 0). Index a runs over the space coordinates. 
There is no summation over a. 

Now we introduce an orthogonal complement to E. This is a set of functions 
fa, (i is 0, x, y, z) which satisfy the conditions: 



(ip, fa) = J jipfadp = 0, 



((pi, (pj) = J -jfafadp = S i:j , (31) 

where we have introduced the scalar product (•, •) to be used below, and where 
5ij is the Kronecker delta, and / is the distribution function. We define a 
projector II which maps j into {E}. The projector is: 



n(j)=j-E(&J)&- ( 32 ) 



The final form of the macroscopic equation depends on the choice of the func- 
tions (pi, specifically, on the conditions of normalization and orthogonality. Let 
us demonstrate a few cases. 

Let us first assume that the functions fa satisfy conditions, (fa, (pj) = 5ij, and 
that the average momentum of the fluid, / fpdp, is equal to zero. Then we 
can take these functions in the form: 



00 = Co/, 

(pa = C a p a f. (33) 
where the constants d should be found from the normalization conditions 
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Clj fd V = Cln=l, C„ = 1/VS, 

Now we calculate the convolutions (0j, j) : 

(M = -D + J ^^-dp = 0, (34) 



*«, j) = ; 1 T Pa ^^dp = < 

<Jmnk B T J ops 

— oo 



-Dn/ y/mnkBT, ii a — (3 
0, if a ^ /3. 



Finally, substituting result (34) into equation (32) we obtain: 

--Pf^ p/( *4 (35) 

Thus, equation (28) is represented as follows: 

df(x,t) , p9/(x,*) /-a$(|r-r7|) = 



<9t m dr dp J dr 

d f df(x,t) 1 \ 
= %(^p- + ^T P/(a; ' t) J- (36) 

This is the Fokker-Planck equation describing a diffusion in the momentum 
space. The short-memory approximation behind the present approach is ap- 
parent in the present example. 

The result (36) can be generalized to the case when the average momentum 
u does not vanish. In this case the basis functions (33) are chosen as follows: 



0o = // Vn, 
= p a - mu a j 

It is easy to check that they satisfy the condition (31). Then equation (30) 
reads, 



d/OM) _ D d ( df(x,t) p-mu \ 

~^^diss ~ V dp~ \^P~ + { ' } ■ {67) 



17 



Thus, in this section we derived the diffusion equation from the iV-particle 
Liouville equation. In contrast to a commonly known Fokker-Planck equa- 
tion, equation (36) is nonlinear because the diffusion coefficient depends on 
the distribution function. Near equilibrium, the diffusion coefficient becomes 
constant, and the result is the usual Fokker-Planck equation (with a mean- 
field extension due to Vlasov's term). Then the usual stochastic interpreta- 
tion in terms of the fluctuation-dissipation theorem applies. We remind the 
reader that the stochastic interpretation of nonlinear equations is, in general, 
a difficult problem [18], and that the result obtained here corresponds to the 
lower-order (r 2 ) of approximation of microscopic trajectories. Higher-order 
corrections will cease to have a stochastic interpretation in terms of the usual 
fluctuation-dissipation theorem. 



4 Equations of hydrodynamics for simple fluid 

The method discussed above enables one to establish in a simple way the form 
of equations of the macroscopic dynamics to various degrees of approximation. 
In this section, the microscopic dynamics is given by the Liouville equation, 
similar to the previous case. However, we take another set of macroscopic vari- 
ables: density, average velocity, and average temperature of the fluid. Under 
this condition the solution to the problem (2) is the local Maxwell distribution. 
For the hydrodynamic equations, the zeroth (quasi-equilibrium) approxima- 
tion is given by Euler's equations of compressible nonviscous fluid. The next 
order approximation are the Navier-Stokes equations which have dissipative 
terms. 

Higher-order approximations to the hydrodynamic equations, when they are 
derived from the Boltzmann kinetic equation (so-called Burnett approxima- 
tion), are subject to various difficulties, in particular, they exhibit an instabil- 
ity of sound waves at sufficiently short wave length (see, e. g. [20] for a recent 
review) . Here we demonstrate how model hydrodynamic equations, including 
post-Navier- Stokes approximations, can be derived on the basis of coarse- 
graining idea, and investigate the linear stability of the obtained equations. 
We will find that the resulting equations are stable. 

Two points need a clarification before we proceed further [5]. First, below we 
consider the simplest Liouville equation for the one-particle distribution, de- 
scribing a free moving particle without interactions. The procedure of coarse- 
graining we use is an implementation of collisions leading to dissipation. If we 
had used the full interacting iV-particle Liouville equation, the result would be 
different, in the first place, in the expression for the local equilibrium pressure. 
Whereas in the present case we have the ideal gas pressure, in the iV-particle 
case the non-ideal gas pressure would arise. 
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Second, and more essential is that, to the order of the Navier-Stokes equations, 
the result of our method is identical to the lowest-order Chapman-Enskog 
method as applied to the Boltzmann equation with a single relaxation time 
model collision integral (the Bhatnagar-Gross-Krook model [21]). However, 
this happens only at this particular order of approximation, because already 
the next, post-Navier-Stokes approximation, is different from the Burnett hy- 
drodynamics as derived from the BGK model (the latter is linearly unstable). 



4-1 Derivation of the Navier-Stokes equations 



Let us assume that reversible microscopic dynamics is given by the one-particle 
Liouville equation, 

df df 

m = (38) 

where / = f(r,v,t) is the one-particle distribution function, and index i 
runs over spatial components {x, y, z}. Subject to appropriate boundary 
conditions which we assume, this equation conserves the Boltzmann entropy 
S = —k B j f In fdvdx. 

We introduce the following hydrodynamic moments as the macroscopic vari- 
ables: M = / fdv, Mi = J Vifdv, M 4 = / v 2 fd\. These variables are related 
to the more conventional density, average velocity and temperature, n, u, T 
as follows: 



3nk B T 2 

M = n, Mi = nui, M4 = h nu , 

m 

n=M , u i = Mo 1 M i , T=——(M 4 -M^ 1 M t M l ). (39) 
The quasi-equilibrium distribution function (local Maxwellian) reads: 



( m \ 3 / 2 (-m{v-u) 2 \ . . 

A = "(£iSr) exp l 2k B r )■ (40) 

Here and below, n, u, and T depend on r and t. 

Based on the microscopic dynamics (38), the set of macroscopic variables (39), 
and the quasi-equilibrium (40), we can derive the equations of the macroscopic 
motion. 

A specific feature of the present example is that the quasi-equilibrium equation 
for the density (the continuity equation), 
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should be excluded out of the further corrections. This rule should be applied 
generally: If a part of the chosen macroscopic variables (momentum flux nu 
here) correspond to fluxes of other macroscopic variables, then the quasi- 
equilibrium equation for the latter is already exact, and has to be exempted 
of corrections. 

The quasi-equilibrium approximation for the rest of the macroscopic variables 
is derived in the usual way. In order to derive the equation for the velocity, we 
substitute the local Maxwellian into the one-particle Liouville equation, and 
act with the operator ii k — J Vk ■ dv on both the sides of the equation (38). 
We have: 



dnuk d nksT dnukUj 
dt dr k m drj 



Similarly, we derive the equation for the energy density, and the complete sys- 
tem of equations of the quasi-equilibrium approximation reads (Euler equa- 
tions): 



On 


dnui 


~dt = 


dr-i 


dnuk 


d nksT dnukUj 


dt 


drk m dr-j 


de 
di~ 


d (hk B T 2 

— — — nui + u nu 

OTi \ m 



(42) 



Now we are going to derive the next order approximation to the macroscopic 
dynamics (first order in the coarse-graining time r). For the velocity equation 
we have: 



p 1 ( [ d 2 fo , \ - d(p nUk 

^ = 2 [J ^dr-dV^ ~ S -dM-ti 

where (f>j are the corresponding right hand sides of the Euler equations (42). In 
order to take derivatives with respect to macroscopic moments {M , Mj, M 4 }, 
we need to rewrite equations (42) in terms of these variables instead of {n, u iy T}. 
After some computation, we obtain: 
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R 



1 d (nk B T 



nu k 



2dr j 



m 



du k du^ _ 2 du n ^ 
drj drk 3 dr n ^ 



(43) 



For the energy we obtain: 



Rr 




j J 



5_d_ ( nk 2 B T dT 
2 dr; \ m 2 dr-j 



(44) 



Thus, we get the system of the Navier-Stokes equation in the following form: 



dn 
~dt 

dnuk 



di 



dnui 
dri 

d nk B T 

drk 1TL 

t d nksT 
2 drj m 



dnukitj 
dr-j 
' duk du 



+ 



dri dr k 




5k B T 2 
nui + u nui 



m 



2 du n 

1 

3 dr n 



5 d (nk 2 B TdT 
m 2 dr; 



T 2dr, 



(45) 



We see that kinetic coefficients (viscosity and heat conductivity) are propor- 
tional to the coarse-graining time r. Note that they are identical with kinetic 
coefficients as derived from the Bhatnagar-Gross-Krook model [21] in the first 
approximation of the Chapman-Enskog method [16] (also, in particular, no 
bulk viscosity) . 



4-2 Post- Navier-Stokes equations 



Now we are going to obtain the second-order approximation to the hydrody- 
namic equations in the framework of the present approach. We will compare 
qualitatively the result with the Burnett approximation. The comparison con- 
cerns stability of the hydrodynamic modes near global equilibrium, which is 
violated for the Burnett approximation. Though the derivation is straightfor- 
ward also in the general, nonlinear case, we shall consider only the linearized 
equations which is appropriate to our purpose here. 



Linearizing the local Maxwell distribution function, we obtain: 
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\2irk B T J \n k B T \2k B T 2 J T 
= | (M + 2M iCi + (Jm 4 - M ) (c 2 - l)}e~ c2 , 



(46) 



where we have introduced dimensionless variables: q = t>j/ 'vt,vt = ^j2k B T /m 
is the thermal velocity, M = 5n/n , Mi = 5ui/vt, M 4 = (3/2)(<5n/n + 
ST /T ). Note that <5n, and <5T determine deviations of these variables from 
their equilibrium values, n , and T . 

The linearized Navier-Stokes equations read: 



<9M dMi 



dt 
dM k 

dt 
dM A 



dri 



ldM A r 
i + _ 

3 dr k 4 dr 



d (dM h dM, 2dM, 



5 dMi 

'2~dT t 



+ T 



J 

5 <9 2 M 4 
2 Sr^ri 



<9r\- 



+ 



dr k 3 <9r 



(47) 



Let us first compute the post-Navier-Stokes correction to the velocity equation. 
In accordance with the equation (10), the first part of this term under linear 
approximation is: 



1_ (dr_d£_ \ _ 1_ ^ ( d<pi d<p* \ 
31^ [df df J ) ^WdMidMp) ~ 

=-\S Ck d^dV n c ^ { (Mo + 2M ^ + - Mo ) ( c2 - 1) } e ~ cVc 

5 d d 2 M A 1 d /3 d 2 M d 2 M 4 \ 5 d d 2 M A 

108 dridr s dr s 6dr k \4dr s dr s dr s dr s ) 108 dr k dr s dr s 

1 g <9 2 M 13 d d 2 M A 

8dr k dr s dr s 108 dr k dr s dr s ' 



The part of equation (10) proportional to the first-order correction is: 




5 d d 2 M A 1 d d 2 M 4 

6 dr k dr s dr s 9 dr k dr s dr s 

Combining together terms (48), and (49), we obtain: 
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(2) _ 1 d d 2 M 89 d d 2 M 4 



Mk 8 dr k dr s dr s 108 dr k dr s dr s 
Similar calculation for the energy equation leads to the following result: 

- / c2 aJ^ c,c ' Ck i (M " + mc ' + ~ Mo ) " §) } e ~"" 3c + 

^25 d d 2 Mi 1 /21 9 <9 2 M; 25 9 <9 2 M; \ 19 9 <9 2 M; 



72 dri dr s dr s 6 \ 4 dri dr s dr s 12 <9r; dr s dr s I 36 dr; dr s dr s 



The term proportional to the first-order corrections gives: 



5 ( d 2 dMi\ 25 ( d 2 dM, } 
+ 



6 \dr s dr s dr^ I 4 \dr s dr s dr 



Thus, we obtain: 



Finally, combining together all the terms, we obtain the following system of 
linearized hydrodynamic equations: 



dM 


dMi 


dt 


dri ' 


dM k 


I dM^ 


dt 


3 dr k 



r_d_ ( dM k + dM 1 _ 2dM n§ \ + 
4 drj \ drj dr k 3 dr n kj J 

r 2 [ 1 _g "9 2 M | 89 g d 2 M, \ 
\ 8 <9r fc dr s dr s 108 <9r fc dr s dr s J ' 
<9M 4 _ 5<9M; 5 <9 2 M 4 2 59 / <9 2 0MA 
<9t 2 9r, 2 dri drj 9 \dr.dr. dri I ' 



Now we are in a position to investigate the dispersion relation of this system. 
Substituting Mj = Mj exp(u;t + i(k, r)) (i = 0, k, 4) into equation (51), we 
reduce the problem to finding the spectrum of the matrix: 
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iky g 



.Ik- 2 - -LP 

4^ 12 x 



12 



12 



ikq. 



12 

-±P - -LP 

4^ 12N 



12 



—ik 



z 



12 



12 

-IP _ -LP 

4^ 12^2 



_ A , 89/£ 

f^x ^3 T 10g 

th, y\3^ 108 

fu 

z \3 108 



ifc x (2 + 59 g ) (2 + 59 g ) ^2(2 + 59 g ) 



2 ft 



This matrix has five eigenvalues. These real parts of these eigenvalues respon- 
sible for the decay rate of the corresponding modes are shown in Fig. 2 as 
functions of the wave vector k. We see that all real parts of all the eigenval- 
ues are non-positive for any wave vector. In other words, this means that the 
present system is linearly stable. For the Burnett hydrodynamics as derived 
from the Boltzmann or from the single relaxation time Bhatnagar-Gross-Krook 
model, it is well known that the decay rate of the acoustic becomes positive 
after some value of the wave vector [6,20] which leads to the instability. While 
the method suggested here is clearly semi-phenomenological (coarse-graining 
time r remains unspecified), the consistency of the expansion with the en- 
tropy requirements, and especially the latter result of the linearly stable post- 
Navier-Stokes correction strongly indicates that it might be more suited to 
establishing models of highly nonequilibrium hydrodynamics. 



4-3 Diffusion in the two- component fluid 



In this example we consider a mixture of the particles of two kinds. We deter- 
mine microscopic equations as two independent one-particle Liouville equa- 
tions. Using our general procedure, we obtain a diffusion behavior of the mix- 
ture on the macroscopic level. 

Microscopic equations of motion for the particles of the first and of the second 
kind are: 



at"" Vl ^' (52) 
df'2 df 2 ( . 

-m=-^- (53) 

We denote all variables related to particles of the first kind with index 1, and 
with the index 2 for the second kind, fx = /i(r, v l5 1) and f 2 = f 2 (r, v 2 , t) are 
one-particle distribution functions. 

In order to describe hydrodynamics of the system, we introduce the following 
macroscopic variables: 
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Mi = pi = mi J fidv, 
M 2 = p2 = m 2 J f 2 dv, 

Mi = (pi + p 2 )ui = mi J v u fidv + m 2 J Vnf-ydv, (54) 
M T = \{n 1 + n 2 )k B T = m, f { ^^hdv + m 2 J ^^/ 2 rfv, 

where pi, p 2 are densities; Ui is the i-th spatial component of the average 
velocity of the mixture; and T is the temperature. 

The quasi-equilibrium distribution functions are: 



^-wfeV^pl-^^^, (55, 

m 2 
2ixk B T(v) 




After the convolution of the system (52)-(53) with the operators p\ — J mi-dv, 
p 2 = jm 2 -dv, pui = fm^-dv, p u2 = jm 2 v 2 -dv, p T1 = J m 1 (v 1 -u) 2 /2-dv, 
PT2 = J m 2 (v 2 — u) 2 /2 ■ dv, and after summations, we obtain the system of 
Euler equations for the binary mixture: 



(57) 
(58) 



dpi 


dpiUi 


dt 


dvi 


dp2 


dp 2 Ui 


dt 


dn 


dpu k 


dnk B T < 


dt 


dr k 


3 dnk B T 


3 duink B T 


2 dt 


2 dn 



dpu k u,i 

*r' (59) 

- nk B T^-. (60) 



Let us now calculate the dissipative correction for the density equations. For 
the equation (57) we obtain: 



R _ 1 / d 2 n 1 k B T d 2 p 1 u i u j \ 
1 2 \ dr 2 dr-idr-j J 

_ ld_ (Mk_d_ f MjPi \ _ MkPi_d_ fM jPl \\ 
2dr k \ p dr-j \ p J p 2 dr-j \ p J ) 

_ 1 d ( M k p 2 d ( Mjp 2 \ \ _ 1 d f pi dnk B T p 1 dnUjU k \ 
2 dr k \ p 2 dr-j \ p ) ) 2 dr k \ p dr k p drj J ' 
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Substituting (54) into (61), and performing similar to the above calculations 
for the equation (57), we arrive at the diffusion equations for the binary mix- 
ture under a simplifying assumption T = const, u = 0: 



dpi _ rk B T j 


[ d 2 n\ 


d 


(? 


dn 


dt' 2 1 


[ dr k dr k 


dr k 




dr k 


dp 2 rk B T j 


f d 2 n 2 


d 


( P2 


dn 


dt' 2| 


[ dr k dr k 


dr k 




dr k 



(62) 
(63) 



Diffusion coefficient rk B T/2m coincides with Einstein diffusion coefficient, 
where r has meaning of the average relaxation time. 



5 Hydrodynamic equation for the fluid with long-range interaction 



In this section we derive equations of hydrodynamics from the nonlinear Vlasov 
equation. This example is also interesting from the methodological point of 
view. Usual methods of reduction of the description are applicable mostly to 
a linear microscopic dynamics [8,7,9]. These methods are based on a formal 
solution to the microscopic equation of motion presented in the exponential 
form. This is not directly possible for nonlinear microscopic models. Since we 
avoid integrating microscopic equations of motion, our approach is immedi- 
ately applicable to nonlinear systems without any modifications. 

Let the microscopic dynamics be given by the Vlasov equation (18). Choosing 
again the usual hydrodynamic fields for the macroscopic variables (39), we 
obtain the system of Euler equations enriched by the mean-field terms: 



dn dnui 

dt dr\ 

dnu k d nk B T dnu k Uj 

dt dr k m dr 



+ nF k , (64) 



j 



de d ( 5k B T 2 \ 

— = -— nui + u nui + 2nUiFi, 

dt dr,i \ m J 

where Fi is the i— th spatial component of the mean-field force given by equa- 
tion (19). 

Following the same route as in section (4.1), we compute term by term the cor- 
rection for (64). The first term in the brackets in equation (9) is proportional 
to: 
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dJ *r- v . iV 92 ^ 1 Vt 9 ( d f° F ) 1 9 ^ 

9/ * 3 ' dridrj 1 dri \dvj 3 ) dvi 



+ F—(v— S ]+F—(—F^\ (65) 
1 dvi \ 3 drj J 1 dvi \dvj 3 J' 

where ^ is given by equation(21). 

In order to calculate the mean-field terms in the velocity equation, we act by 
the operator / Vk • dv on equation (65). We obtain: 



-nm k + -Q^-Fi. (66) 



The second part of equation (9) is: 



dcf) u . , ,. , , i)M, n dFi ^ dM, Mi <9F 7 M „ dM, 
^ dMj drj drj drj M drj drj 

where Mj correspond to (39). Rewriting this expression in terms of the vari- 
ables n, u and T, and combining the result with equation (66), we obtain: 



R { nl k = n* k . 



Now let us calculate the correction to the energy equation (45) in the presence 
of the mean-field interaction. Action of the operator fv 2 -dv on equation (65) 
gives: 



d ( k B Tn\ dn , dn \ 

~ 5 TT~ ~~ 2u i u 3T> u TT~ \ F 3 ~ 

v drj \ m J or, i orj J 

' duj dui dui\ 

2nui—^- - 2nui- 2nu j — Fj + 

v ori drj dr\ J 

' ^ nksT x 9 2 A \ dFj 

—5 da — ZnUiUj — u ndij — J - 

m dn 



d fk B Tn\ dnu 2 ^ 
' dr-i \ m I dr-i 



-2n Ui % + 3— — + — — Fi + 2nF t F l . 



The differential term for the energy density equation gives: 
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^ dcf> £ M 4 8M dF t 2 9// M n M, 



dM; M <9ri "<9r; 3 dn \ V M 

2FjMj- — — * - 2F i M i — — * - 2- J 



4 J 9rj VM 7 ~ J ~'~ l drj \M J ~ M dr,- 
9 / M n M„\ Fj / M n M„\ <9M 

" ^ I 4 ~ ~m^) + Mo I 4 " "l^T J ~d7~ 

2M j F j dM, _ p 2 8_ / _ M n M n 



M dn l 3dn V M ( 



'0 



.7 



g /MA 2FjMj dM. 
1 j d rj \Mj Mo dr 



Thus, we obtain the first-order correction to the energy equation, 



i# } = — —Fi— + 2F i u i u j — - 2nui*i + 2nF^, 
m or,- or., 



Finally, we arrive at the following system of hydrodynamic equations with the 
the mean-field interaction: 



<9n dnuj 

i 

h.r/V Bn.1l.i-1l.: 

+ nF k + 



dt dn 
dnuk d nksT dnu^Uj 



dt drk m dr-j 



t d nk B T ( du k duj 2du n \ rn 




^ = - ( —^-nui + w 2 nw ; | : 2iin t h) 
dt dri \ m 

5 d (nk 2 B TdT\ 3k B T^dn „ dn 

T o7T 2~^~ + T Fi TT + rFiUiUjT. 

2 dri \ m A dri J m dr { drj 

rnu^i + rnFiFi, 



Equations (67) are the general form of the hydrodynamic equations of a simple 
fluid with the mean-field interaction. Examples of the system for which this 
result may be relevant can be found in studies of electron transport in the 
various media [19], as well as in description of non-Newtonian fluids [22]. For 
each particular case the interaction potential $(|r — r'|) has its specific form, 
and leads to the corresponding hydrodynamics of the system. 
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6 Conclusion 



In this paper the formalization of Ehrenfest's approach to irreversible dynam- 
ics is given in details. This method allows one to derive macroscopic equations 
of motion on the basis of the microscopic dynamics and the very transparent 
coarse-graining procedure. The method is applicable to both reversible as well 
as to the irreversible microscopic dynamics, independently of whether it is 
linear or not. We have presented a set of examples demonstrating how this 
method is applied to various situations. 

The most interesting continuation of this approach is, of course, how to specify 
in a sensible and practical way the coarse-graining time r in order to make 
the modelling parameter-free. This requires is a subject of our current studies 
(see [23,24,25]). 

Finally, whereas we have focused on the application of our formalism to the 
entropy- conserving microscopic dynamics, it should be mentioned that it is 
applicable also to constructing slow invariant manifolds of dissipative systems. 
In particular, when applied to the Boltzmann equation, the result is equivalent 
at r — > oo to the exact Chapman-Enskog solution [23]. 
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Fig. 2. Attenuation rates of various modes of the post-Navier-Stokes equations as 
functions of the wave vector. Attenuation rate of the twice degenerated shear mode 
is curve 1. Attenuation rate of the two sound modes is curve 2. Attenuation rate of 
the diffusion mode is curve 3. 
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